Draft version January 12, 2013 

Preprint typeset using I^-T^X style emulatcapj v. 04/21/05 
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Abstract 

The recently-discovered accreting X-ray pulsar IGR J17480-2446 spins at a frequency of ~11 Hz. We 
show that Type I X-ray bursts from this source display oscillations at the same frequency as the stellar 
spin. IGR J17480-2446 is the first secure case of a slowly rotating neutron star which shows Type I 
burst oscillations, all other sources featuring such oscillations spin at hundreds of Hertz. This means 
that we can test burst oscillation models in a completely different regime. We explore the origin of 
Type I burst oscillations in IGR J 17480-2446 and conclude that they are not caused by global modes 
in the neutron star ocean. We also show that the Coriolis force is not able to confine an oscillation- 
producing hot-spot on the stellar surface. The most likely scenario is that the burst oscillations are 
produced by a hot-spot confined by hydromagnctic stresses. 

Subject headings: pulsars: individual (IGR J17480-2446) — stars: magnetic field 
— stars: rotation — X-rays: bursts 



stars: neutron 



1. INTRODUCTION 

Accreting neutron stars (NSs) in low mass X-ray bi- 
naries show bright Type I X-ray bursts. These begin 
with a rapid increase of the X-ray flux (the rise) followed 
by a slow decrease (the tail) to the pre-burst luminos- 
ity, and last ~10-100 s. These bursts are powered by 
thermonuclear runaways, which burn up a layer of accu- 
mu lated light elements on the N S surface (for a review 
see IStrohmaver fc Bildstenl I2003D . 

A significant fraction of bursts display quasi- 
periodic modulati o ns, kn own as burst oscillations (BOs, 
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I1996T ) . In the discovery paper, 
suggested that BOs were related to the 
spin frequency of the NS. This has been confirmed in five 
accretion-powered millisecond X-ray pulsars (AMXPs), 
where it was found that BO freq uencies are within a few 
perce nt of the spin frequencies (IStrohmaver fc Bildstenl 
2003!). This implies that BOs are caused by a near- 
stationary temperature asymmetry which persists in the 
surface layers of the star during burst. 

The detailed phenomenology of BOs, however, is di- 
verse. In the sources which have evidence for substan- 
tial magnetic fields (the persistent AMXPs), for example, 
BOs appear at a near constant frequency in the tail, with 
some fast chirps in the rise (ICh akrabart v et all 120031: 
IStrohmaver et al.ll200l lAltamirano et al.ll2010al) . In the 
intermittent AMXPs, and NSs without evidence for a dy- 
namically important magnetic field, BOs t ypically drift 
upwards by a few Hertz during the burst (jMuno et al.l 
l2002at I Watts et al.ll2009D . Both the origin of the sur- 
face temperature asymmetry that causes the BOs, and 
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the reason for the observed frequency drifts, remain un- 
solved puzzles. 

One possibility is that the asymmetry is caused by 
global modes (waves) that develop in the bursting 
ocean ([Strohmaver fc Led fl996t iHevH l200i iCumrnind 
I2005T ). As the ocean cools, its scale height H de- 
creases and the pattern speed, which scales as \H, 
changes, leading to frequency drift. The ki nd of 
mode s that mi ght be excited include r- modes (jHevll 
120041) , g-modes (|Bildsten fc Cumming 1998) or magncto- 
hvdrodvnamical modes ([Heng fc Spitkovskvl 120091) . To 
date none of these models have managed to explain both 
the obse rved frequencies and the magnitude of the drift s 
(see e.g. iPiro fc Bildstenl l2005t iBerkhout fc Levinll2008| ) . 

An alternative possibility is that a compact burn- 
ing h ot-spot develops on the surface ([Strohmaver et all 
119961 ). The question is then how confinement, of fuel 
or the flame front itself, might be achieved. For un- 
magnetized stars, the Coriolis force could be a n effective 
confining mechanism (fSpitkovs kv et all 120021 hereafter 
SLU02). Although lifting and expansion of the hot fluid 
should cause spreading of the burning fuel, the Coriolis 
force would oppose such motion of the flame front by de- 
flecting its velocity. This mechanism could account for 
the presence of oscillations in the burst rise, although it 
does not easily explain the presence of BOs in the tail 
or the frequency drifts. The Coriolis force was however 
attractive in that it could explain why BOs had not been 
seen in any NS with spin frequency v s < 245 Hz (since 
for more slowly rotating stars the Coriolis force is not 
dynamically relevant, see Sec. 14. 2j) . 

A strong magnetic field could also lead to confinement 
of fuel or flame, with the restoring force supplied by field 
pressure or stress. This mechanism is particularly plau- 
sible for accreting pulsars, where the existence of dy- 
namically important magnetic fields is suggested by the 
presence of accretion-powered pulsations ( APPs) . 

Since there is no model that definitively explains all 
features of BO phenomenology, it is important to explore 
the applicability of the models under the widest possi- 



2 



ble range of conditions. Until now, however, the only 
stars to show BOs were rapid rotators with v s > 245 Hz. 
This situati on has changed with our discovery of burst 
oscillations (|Altamirano et a l. 2010b) from the accreting 
pulsar IGR J17480-2446 (hereafter J17480), which ro- 
tates an order of magnitude m ore slowly at ~ 11 Hz 
(jStrohmaver fc Markwardtll2010D . 

2. OBSERVATIONS AND DATA ANALYSIS 

J17480 was detected in outburst on 10 October 2010 
in a Galactic bulge scan with INTEGRAL (jBordas et al.1 
[20101) . The source is located in the globular cluster 
Terzan 5 and its outburst lasted for ~55 days (ac- 
cording to MAXI monitoring observations) before the 
source became undetectable due to solar constraints. 
Further monitoring observations by MAXI showed that 
the source had probably returned to quiescence by Jan- 
uary 2011. J17480 is an 11 Hz pulsar in a 21.3 
hr orbit around a companion with M > 0.4 
(jStrohmaver fc Markwardtl 120101: iPapitto et all 1201 ll ). 
Follow-up observations of RXTE /PC A have collected 
approximately 294 ks of data to 19 Novemb er 2010. We 
used all RXTE /VGA (jJahoda et all 120061 ) public ob- 
servations of J17480 taken between MJD 55482.0 and 
55519.2 (October 13-November 19, 2010). For our co- 
herent timing analysis we used Event and GoodXenon 
data modes, keeping only absolute channels 5-37 (~ 
2 — 16 kcV), rcbinned to 1/8192 s and bary centered with 
the FT OOL f axbary using the position of iPoolev et al.l 
< l20Toh . 

2.1. Lightcurve and accretion rate 

We estimate the flux during each observation from 
the background subtracted 2 —16 keV Crab normalized 
intensity (estimated following Ivan Straaten et"alll2003l) 
and converted it to erg cm -2 s _1 assuming a Crab spec- 
trum for the source. Flux is converted into lumi nosity 
using a distance of 5.5 kpc (jOrtolani et al.l [20071). Lu- 
minosity L x is converted into mass accretion rate M 
using L x = 20%Mc 2 ([Frank et al.ll2002D . The value 
of M obtained is approximate, since the X-ray flux 
is not bolomctric, and no correction has been applied 
for absorption, the unknown disc inclination or gravita- 
tional redshift. Conversion of bolometric luminosity to 
mass accretion rate is further complicated by the un- 
certainty in the radiative efficiency. The overall uncer- 
tainty on M is within a factor ~ 3, taking into account 
all these corrections (jin't Zand et al.ll2007f ). The low- 
est luminosity observed gives M = 9 x 10~ 10 M Q yr _1 
(in the first observation) while at the outburst peak 

M = 7 x 10" 9 M Q yr" 1 , i.e., 5% and 37% of the Edding- 
ton rate MEdd = 1-9 x lO _8 M0yr _1 (assuming solar 
composition for the accreted gas). 

2.2. Bursts 

Bursts were identified visually in a 1 s lightcurve. The 
beginning and end of the bursts were defined as the 
points where the flux first and last exceeded the max- 
imum pre-burst flux in the same observation. Typ- 
ical burst length is ~ 100 — 200 s. Burst recur- 
rence time and burst peak-to-persistent flux ratio de- 
crease with increasing persistent flux until the out- 
burst peak when the bursts disappear and millihcrtz 
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Fig. 1.— Top panel: 2-16 kcV X-ray lightcurve of J17480 av- 
eraged over 16 s long data intervals (red dots) and average flux 
per Obsld (black crosses). S haded regions contain the bursts with 
clear evidence of cooling (Chakraborty & Bhattacharyya] 120111 ; 
ILinares et al J 1201 lft . note that the first region has a different scal- 
ing for both axes. Vertical dashed lines mark the region where 
X-ray bursts are difficult to define, so are excluded from the burst 
analysis. Middle panel: Sinusoidal fractional amplitudes of the 
APPs for the fundamental frequency and first two overtones. Frac- 
tional amplitude is ~ 28% in the first observation, subsequently 
dropping to 1 — 2%. The fundamental and first overtone are de- 
tected throughout the outburst; overtones are seen only sporadi- 
cally. Bottom panel: Sinusoidal fractional amplitude of the BOs 
for the fundamental frequency and first two overtones. Fractional 
amplitude evolution is very similar to that of the APPs, dropping 
after the first observation from 28% to ~ 1 — 4%. 

quasi-periodic oscillations ([Linares et al.l I2010D appear. 
As the outburst flux decreases again, the bursts grad- 
ually reappear with recurrence time and burst peak- 
to-persistent flux ratio increasing wi th decreasing flux 
([Chakraborty fc Bhattacharvval 1201 ll) . Between MJD 
55486 and MJD 55492 the difference between bursts and 
persistent flux fluctuations becomes negligible and burst 
recurrence time is of the same order as the duration of the 
bursts themselves, therefore burst rise and tail become 
difficult to define. In this case we do not identify the 
bursts, and for the purpose of measuring the oscillations 
treat all data as persistent flux. 

Using this identification criterion, we found 231 bursts. 
As discussed by other a uthors, some showed clear spec- 
tral evidence for cooling (TChakrabortv &: Bhattacharvval 
120111: iMotta et all 120111: ILinares et all 120111) . were ob- 
served when the persistent luminosity was < 15% Ed- 
dington and were identified conclusively as Type I (ther- 
monuclear) bursts (Fig. [TJ. The nature of the other 
bursts, observed when the persistent flux was ~ 15 — 
35% Eddington, remains debated due to the lack of 
clear evidence for cooling (iGallowav fc in't Zand! 120101 : 
iChakrabortv fc Bhattacharvvall201lD . 

2.3. Pulsation analysis 

The data were folded in non-overlapping intervals of 
~ 10 — 20 s for the first observation (where APP and 
BO fractional amplitudes are ~ 30%), and ~ 100 — 500 s 



3 



1 1 .5 1 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 1 1 200 




200 



10.0 

50 100 150 200 

TIME (S) 

Fig. 2. — Light curve and dynamical power spectrum for the 2010 
October 13 thermo nuclear (as showed by the spectral analysis of i.e. 
ILinares ct al. 20lJ) burst, using GoodXenon data from PCU2, the 
only detector active at the time. The dynamical power spectrum 
uses overlapping 4 s bins, with new bins starting at 1 s intervals. 
The contours show Leahy normalized powers of 30, 60, 90, 120 and 
150. 



TABLE 1 





v [Hz] 


!> [10~ i2 Hzs- 1 ] 


Fundamental BOs 
Fundamental APPs 


11.04488532(3) 
11.04488540(5) 


1.44(3) 
1.42(5) 


1st Overtone BOs 
1st Overtone APPs 


11.04488547(5) 
11.04488548(10) 


1.33(5) 
1.34(11) 


2nd Overtone BOs 






2nd Overtone APPs 


11.04488530(8) 


1.40(6) 



Note. — All frequencies, inferred from each harmonic, refer 
to the epoch MJD 55482. In brackets are the statistical errors 
at the 68% confi dence level, calculated by means of Monte Carlo 
simulations as in IHartman etaLI pOOM ). 

for all others using the ephemeris of lPapitto et al.l (|2011l ) 
to search for both APPs and BOs. This means that in 
some cases (when the signal is sufficiently strong) mul- 
tiple BO profiles are constructed for each burst. Only 
pulse harmonics with SNR> 3.5 were retained (i.e., har- 
monics whose amplitude was larger than 3.5 times the 
amplitude standard deviation), giving less than 1 false 
pulse detection. 

Each pulse profile was decomposed into the funda- 
mental frequency and overtones to compute fractional 
amplitudes and phases (we were able to identify up to 
five overtones in some intervals). The pulse phases of 
each profile were calculated for each harmonic and an- 
alyzed separately. T he method applied is described in 
iPatruno et al.l (|2010D . along with the procedure to cal- 
culate pulse amplitudes and statistical errors. Phases 
were then fitted with a Keplerian orbit, a linear and a 
parabolic term representing the pulse frequency and its 
first time derivative. 

3. DISCOVERY AND PROPERTIES OF BURST 
OSCILLATIONS 

Burst oscillations were detected in all 231 bursts iden- 
tified. Fig. [T] shows the fractional amplitudes of the 
BOs compared to those of the APPs. The fractional am- 
plitudes of BOs are always of the same order or larger 
than those of APPs, and always significantly larger than 
expected if the modulations stem only from residual per- 
sistent emission. The detection of BOs can therefore be 
considered secure. 



Note that the APP contribution has not been sub- 
tracted when calculating BO fractional amplitudes and 
phases. The measured BO fractional amplitudes are 
comparable to or hig her than the APP ones, it can be 
seen from E q (11) of IWatts et all ((20051 ) and Eq. (1) of 
iWatts et~a l. (2008) that the correction is less than a few 
percent. 

The presence of tim ing noise, as seen in other AMXPs 
(|Patruno et al.l I2009D complicates timing analysis of 
APPs and BOs. The scatter observed in the oscillation 
phases has an amplitude of 0.1-0.2 cycles within each 
Obsld for both BOs and APPs. However adjacent BOs 
and APPs are always phase coincident and phase locked 
to within ~ 0.05 cycles. In the first observation, for ex- 
ample, which has the highest SNR, the phase difference 
between the peak of the fundamental for APPs and BOs 
is 0.014 ± 0.005 cycles (5.1 ± 1.8 deg). This coincidence 
of results is remarkable given the strong timing noise, 
highlighting the fact that BOs do indeed closely track 
APPs. 

To check whether BO phases have the same temporal 
dependence as the APPs, we fit a pulse frequency model 
plus its first derivative to the the entire two data sets of 
BO and APP phases separately. The pulse frequency and 
derivative obtained are consistent with being the same to 
within two standard deviations for both the fundamental 
frequency and the first overtone (Tab. [T]). 

Inspecting the behavior of the pulse frequency locally 
(comparing BO and APP pulse frequencies in each indi- 
vidual Obsld), pulse frequency is again consistent with 
being the same for the two data sets within to two 
standard deviations. In the first observations, when 
the BO and APP amplitudes are around 30% (Fig. [2), 
the frequencies of APPs (11.044881(2) Hz) and BOs 
(11.04493(7) Hz) are identical within the errors (~ 7 x 
10" 5 Hz). 

Comparing the properties of J174 80 to those of the 
other persistent pulsars with BOs (IChakrabarty et al.l 
2003t IStrohmaver et al.l 120031: lAltamirano et al.l l2010aF 
Riggio et al.M201lD . JT7480 shows BOs in every burst, 
just like SAX J1808.4-3658, XTE J1814-338 and IGR 
J17511-3057. J17480 is most similar to J1814: BOs are 
present throughout the bursts, BOs and APPs are phase- 
locked very closely, BOs have strong harmonic content 
and no measurable frequency drifts. J1808 and J17511, 
by contrast, have BOs with weaker harmonic content, 
which do not persist throughout the bursts, and display 
fast drifts in the rise. In terms of BO amplitude, the 
persistent pulsars differ: J1808 has amplitudes of 3-5% 
r.m.s., J17511 of 5-15% and J1814 of 9-15%. J17480 has 
BO amplitude of 30% r.m.s. in the first burst, dropping 
to ~ 2% thereafter. However, in all cases the BO am- 
plitude is comparable to (within a few percent) the APP 
amplitude. 

4. BURST OSCILLATION MODELS FOR J17480 

We first show how we can exclude mode and Coriolis 
force confinement models based on the coincidence of BO 
and APP frequencies and on slow rotation. 

4.1. Global Modes 

Given a mode with azimuthal number m and frequency 
u r in the rotating frame of the star, the frequency v a an 
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inertial observer would measure is 

v = mv s + v T (1) 
the sign of v x being positive or negative, depending on 
whether the mode is prograde or retrograde. 

In the first, and best constrained burst (Section [3]), v 
and v a were found to differ by no more than 10 -4 Hz. 
Excluding modes with m > 2, which cannot explain the 
high f ractional am plitudes (~ 30%) observed in the first 
burst ()Hevll[2004T ). leaves us with two main possibilities. 

The first is that we have modes with |z/ r | ~ v s and 
m=0 (prograde) or m=2 (retrograde). Even though 
modes with frequencies arou nd 10 Hz, such as g-modcs 
(|Bildsten fe Cumminel [1998) do exist, it would require 
extreme fine tuning of parameters for v T to match v a 
within ~ 10 -4 Hz. We are not aware of any mecha- 
nism that could tie the frequency of a mode so closely to 
the spin. 

The second option would be to have prograde modes 
with m = 1 and v x ~ 10~ 4 Hz. However the process that 
excites a mode must have a time scale r ~ \ jv. In this 
case r ~ 10 4 s, which is longer than the duration of the 
bursts. We therefore conclude that global modes cannot 
explain the presence of BOs in J17480. 

4.2. Coriolis force confinement 

For rapidly rotating stars. ISLU02I showed that burning 
could be confined, due to the effect of the Coriolis force, 
over a length sca le of order the Rossby adjustment radius 
(lPedloskvlll987h : 



R R = ^H/Anv a = 34 km m\'2 R w ' H^ 2 u^ (2) 
where Mi. 4 is mass in units of 1.4 M Q , Rio radius in 
units of 10 km, Hip the s cale height of the burning fluid 
in units of 10 m (|SLU02f ) and u e>10 = v s /{10 Hz). 

For J17480, i?R significantly exceeds the stellar radius: 
Coriolis force confinement is not effective. So although 
this mechanism might be important for faster NSs (for 
v 8 > 200 Hz, i?R < 1.7 km), it cannot cause BOs in 
J17480. 

5. MAGNETIC CONFINEMENT 

Having ruled out these two models, we now explore 
the possibility that the magnetic field could confine the 
hot-spot. We begin by placing limits on the strength of 
the magnetic field for this source, following the stan dard 
procedure outlined in lPsaltis fc Chakrabartvl (|1999f ). For 
the magnetic field to be strong enough to channel the 
accretion flow, and hence generate APPs, the magne- 
tospheric radius (the point at which the magnetic field 
disrupts the disc) must exceed the stellar radius. This 
gives a lower limit on the magnetic field. In addition, to 
avoid propeller effects (which would inhibit accretion), 
the magnctosphcric radius has to be lower than the coro- 
tation radius. This gives an upper limit on the magnetic 
field. Using the accretion rates derived in Section 12.11 
we conclude that the magnetic fi eld for J 17480 is in the 
range 2 x 1 8 to 3 x 10 10 G fseelPsaltis fc Chakrabartvl 
(|1999h and lAndersson et al.l ()2005f ) for a detailed expla- 
nation of how to calculate the magnetic field). 

We now need to estimate whether this field is strong 
enough to confine the burning material in order to gen- 
erate BOs. We start by considering static fuel confine- 
ment. A mountain of accrete d material induces a pres - 
sure gradient within the fuel. iBrown fc Bildstenl (|1998l ) 



computed the field that can statically compensate for 
such a gradient by magnetic tension. Confinement within 
an area of radius R Cl at the ignition column-depth ~ 10 8 
g/cm 2 requires magnetic field 



B>3x 10 1U G 



R, 



1 km 



-1/2 



(3) 



Fuel confinement may however be ineffective due 
to MHD intercha nge or "ballooning" instabilities 
(|Litwin et al.1 120011 ) , which allow accreted fluid to es- 
cape the polar cap region and spread over the NS sur- 
face. Nonetheless, confinement might still be achieved 
dynamically by the motion of the fluid itself. This is 
because such motion can induce a horizontal field com- 
ponent which can act to halt further spreading. This 
mechanism requires a weaker field 



B > 4 x 10 9 G 



1 km 



(4) 



dHeng fc Spitkovskvl[2009l) . 

Both of these estimates are compatible with the mag- 
netic field inferred for this source from its accretion 
properties. We conclude that magnetic confinement of 
the burning fluid is a viable mechanism to explain the 
BOs in J17480. Detailed theoretical studies will how- 
ever be needed to verify this possibility, since it depends 
on poorly understood details of the interaction between 
the magnetic field and the hydrodynamics of the burning 
ocean. That said, confinement of burning material at the 
polar cap provides a natural explanation for the phase- 
locking and coincide nce of the APPs and BOs (see also 
lLovelace et al.|[2007h and would explain why the ampli- 
tudes of the two sets of pulsations track each other as 
the accretion footprint varies during the outburst (Fig. 

6. CONCLUSIONS AND WIDER IMPLICATIONS 

We have analyzed the bursts of the slowly rotating X- 
ray pulsar J17480 and found oscillations at the same fre- 
quency as the NS spin. We have shown that neither 
global modes nor a Coriolis force confined hot-spot can 
explain the presence and frequency of these BOs. We 
suggest that the magnetic field could potentially provide 
the necessary confining force for a hot-spot to persist. 
The field needs to be at least B > 10 9 G and this re- 
quirement is compatible with the constraints set by the 
accretion process. This model would neatly explain the 
phase-locking of APPs and BOs, and the fact that the 
amplitudes of the two pulsations track each other. 

We noticed that among the other persistent AMXPs 
XTE J1814-338 is the most similar one, showing almost 
constant frequency of the BOs, harmonic co ntent and 
remarkably phase locking. IWatts et al.l (|2008h excluded 
the possibili ty of magnetic con fi neme nt, since using the 
estimates of lBrown fc Bildstenl (|1998t) (Eq. the field 
was not strong enough to confine the fuel. However, they 
did not consider the case of dynamical confinement when 
spreading is allowed (Eq. H|) . This operates at lower fields 
which are compatible with the magnetic field estimated 
for J1814. 

This mechanism seem less likely for other AMXPs 
such as SAX J1808.4-3658 and ICR J17511-3057, 
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where the fields are lower (|Hartman et al.l I2009T) 
and there are frequency drifts in the rising phase 
of bu rsts (|Chakrabartv et aT] 120031 : lAltamirano et al.l 
l2010a|) . That said, magnetic confinement would offer 
an explanation for the presence of higher harmonics in 
pulsar BOs, since the emitting are a would be bounde d 
(probably near the rotational poles. [Muno et ai1l2002bl ). 

What about the non pulsating sources and the inter- 
mittent pulsars with a weaker magnetic field? They do 
not show harmonic content and their BOs show drifts 
during the decay of the bursts. If magnetic confine- 
ment were responsible one would expect to see APPs 
from these sources as well. Unless the field can be 
boosted temporarily to detectable levels only during the 



burst (jBoutloukos et all I2010D. a different mechanism 
may be required ([Watts fc Strohmaverl [200 6^ . If global 
modes or Coriolis confinement are responsible, then the 
wider emitting surface involved, or the involvement of 
higher colatitudes, would suppress t he presence of har- 
monics, in line with the observations ([Muno et~aT1l2002bl : 
IStrohmaver et al.ll2003l ). 
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